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We study the deformation of a Fermi surface (FS) near charge-ordering (CO) transition. By applying a 
fluctuation-exchange approximation to the two-dimensional extended Hubbard model, we show that the 
FS is largely modified by strong charge fluctuations when the wave number of the CO pattern does not 
match the nesting vector of the FS in a noninteracting system. We also discuss the enhanced anisotropy in 
quasiparticle properties in the resultant metallic state. 

KEYWORDS: organic conductors, charge ordering, Fermi surface, anisotropic electron transport, magne- 
toresistance 
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The shape of a Fermi surface (FS) is an important factor in 
determining the electronic properties of a metal. For many or- 
ganic conductors, band structures and FSs have been obtained 
by the extended Hiickel method 1 ' 2 ' and first-principles calcu- 
lations. 3 ' 4) The band structures thus obtained provide a basis 
for constructing effective tight-binding models, which explain 
the various electronic phases in different types of organic 
salts in a unified way. 51 The shape of the FS is probed with 
high sensitivity by transport measurements such as angle- 
dependent magnetoresistance oscillations, the Shubnikov-de 
Haas effect and magnetoresistance. 6 ' 

However, the shape of the FS is modified by strong elec- 
tronic exchange-correlation effects, which are neglected in 
band calculations. Deformation of the FS has already been 
discussed on the basis of a single-band Hubbard model to 
study high-r c superconductors. Theoretical studies using the 
second-order perturbation theory, 7 ~ 9 ' the one-loop approxi- 
mation 10 ' and the fluctuation-exchange (FLEX) approxima- 
tion 11 ' 12 ' indicate that the FS is gradually deformed as antifer- 
romagnetic (AF) spin fluctuations develop near the AF tran- 
sition. In these studies, the FS is deformed so that its nesting 
condition improves at the wavenumber of the AF ordering. 
However, the degree of deformation obtained is small. 7 ~ 12 ' 

In organic conductors with a 3/4-filled band, charge fluc- 
tuations develop near the charge-ordering (CO) transition and 
strongly influence the electronic properties in a manner differ- 
ent from that of AF spin fluctuations. 13 ' The effect of charge 
fluctuations has been clarified in recent theoretical studies of, 
for example, superconductivity mediated by charge fluctu- 
ations, 14_17 ' non-Fermi-liquid behavior, 18,19 ' anomalous en- 
hancement of Pauli paramagnetism 20 ' and instability toward 
inhomogeneous electronic states. 21 ' Although the electron 
exchange-correlation effect is also expected to modify the FS 
shape via charge fluctuations, such an effect has not yet been 
studied. 

In this paper, we study how charge fluctuations affect the 
FS shape. By employing a FLEX approximation, we demon- 
strate that in contrast to AF spin fluctuations, FS is largely 
modified near the CO transition. We emphasize that this re- 
markable change in FS originates from the large discrepancy 
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Fig. 1 . (a) Schematic of the extended Hubbard model with hopping integral 
t and intersite Coulomb interactions V and V considered in this paper and 
(b) its relation to the crystal structure of 8-ET salts. 



between the CO wave vector and the nesting vector in a nonin- 
teracting system. As a result of FS deformation, quasiparticle 
properties become more anisotropic near the CO transition. 
We also discuss the experimental relevance of our results. 

We consider the extended Hubbard model on a two- 
dimensional square lattice. The Hamiltonian is 
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where c. (c lcr ) is the creation (annihilation) operator of an 
electron on site i with spin <x = f or X, « (cr = c! tv and n,- = 
ii/i + n/i. Here, t is the hopping integral between neighbour- 
ing sites, and U and V are the onsite and nearest-neighbour 
Coulomb interactions (denoted with (i, j)), respectively. In the 
last term of Eq. (1), we have added an intersite Coulomb re- 
pulsion V between next-nearest-neighbour sites in one diag- 
onal direction, as shown in Fig. 1 (a). 

The present Hamiltonian (1) is an effective model for CO 
in the series 0-ET 2 MM'(SCN) 4 [ET=BEDT-TTF, M=T1, Rb, 
M'=Co, Zn] (abbreviated as 0-ET salts), 22 ' 23 ' which exhibits 
the quasi-two-dimensional molecular arrangement shown in 
Fig. 1 (b). Owing to the crystal structure of 6>-ET salts, the 
overlap between molecular orbits is approximately restricted 
to the x and y directions, whereas strong intersite Coulomb 
repulsion remains along one diagonal direction. The Hamil- 
tonian (1) has been adopted as a minimum model in several 
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Fig. 2. Static charge susceptibility xAl^ 0) for U = 5 and kgT = 0.1. The 
CO transition occurs at V = 2.523. Q = (2n/3, 2tt/3) is the wave vector 
characterizing the CO pattern. 



theoretical studies of the CO phenomena in f9-ET salts. 24 30) 

We introduce here a FLEX approximation 31 ' to study the 
FS deformation caused by intersite Coulomb interactions. The 
one-particle Green's function G(k) is related to the self-energy 
Z(fc) through the Dyson equation 

GW = ■ —~ 1 ™y (2) 

;w„ + p — en — Zu(k) 

with a combined notation £ of the wave number k and the 
fermionic Matsubara frequency iw„. Here the noninteracting 
band dispersion is given by = 2?(cos k x + cos k y ) and the 
Hartree term is absorbed in the chemical potential as p. = fx - 
(U/2 + 4V + 2V')n with the average number;? of electrons per 
site. In the FLEX approximation, the self-energy is expressed 
as 



G(k - q) [V c (q) + V s (q)] ■ 
s(q) 



1 + v c 



(3) 



(4) 



Aq)x(q) 

where V c , s {q) represents the exchange-correlation interaction 
potential, v c {q) — U + 4V(cosg A + cosq y ) + 4V cos(q x + 
q y \ v,(f) = -U, X (q) = ~ J k G(k + q)G(k) and _j = 

T Tioj„ J_ J_ &qlQ-n) 2 give the combination of the sum and 
integral with respect to the bosonic Matsubara frequency and 
momentum. We omitted the particle-hole and particle-particle 
ladder diagrams, because SU(2) symmetry of the spin sector 
of electrons is not relevant if the system is far from the AF in- 
stability. 15 - I6) This simple approximation is sufficient to study 
nontrivial exchange-correlation effects on the FS. 

Throughout this paper, we fix the hopping integral, the on- 
site interaction, and the temperature as t = 1, U = 5, and 
k%T = 0.1, respectively, and choose V = V for simplicity. 
We checked that no AF instability occurs in the entire param- 
eter region considered. By combining Eqs. (2)-(4), we self- 
consistently determine G(k) at 3/4 filling (n = 3/2) with a 
relative precision of 10~ 6 . We adopt the units e = % = 1, 

We begin by showing in Fig. 2 the static charge suscepti- 
bility 

X(q,0) 




Xc(q,0) = 



i+v c (?M?,0)' 



(5) 



Fig. 3. (a) V dependence of FS. One quarter of the FS is plotted in the main 
panel with the entire FS shown in the inset, (b) Contour plot of Re I. R (k, 0) 
for V = 2.43. 



As V(= V) increases, the charge susceptibility increases at 
a specific wave vector Q = (2^/3,2^/3), and diverges at 
V = V C = 2.523. The divergence of Xc(Q> 0) indicates a phase 
transition into three-fold-type CO. 26 ' In the following discus- 
sion, we focus on the region V < V c , in which charge fluctua- 
tions become large enough to modify the shape of the FS. 
In Fig. 3 (a), we show a FS determined by 

s k „ +ReI s (Jfc F ,0)-/i = 0, (6) 

where kf is the Fermi wave vector and l, R (kf, 0) is obtained 
from the analytic continuation of 2.(kf, io)k) in the upper plane 
(wk > 0) by using the Pade approximant. For weak inter- 
site Coulomb repulsion (V = 1), the FS (dashed line) is 
nearly circular, and is almost the same shape as that of the 
noninteracting electron system (dotted line). As V(= V) in- 
creases, the FS shape gradually changes. Near the CO tran- 
sition (V = 2.43), the FS (solid line) shape clearly departs 
from that of the noninteracting system. The wave vector Q of 
the CO pattern is also shown by the thick arrow in the inset 
of Fig. 3 (a). The FS shape is clearly modified to assist the 
nesting condition; a fiat part of the FS, which is spanned with 
the wave vector Q, grows as the system approaches the CO 
transition. This is the main result of this study. 

In the present study, the CO wave vector Q is a large mis- 
match for the nesting vector, which is determined by the shape 
of the FS for a noninteracting system, because Q is deter- 
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mined not only by ^(q) but also by v c (q) — U + 4V(cos q x + 
cos q y ) + 4V cos(q x + q y ). However, because Qaf is not af- 
fected by the bare spin interaction v s {q) — —U and is deter- 
mined only byxiq, 0) [Eq. (4)], the wave vector Qaf of the AF 
fluctuations already properly spans the FS of the noninteract- 
ing system. As a result, the modification of the FS induced by 
AF spin fluctuations is small compared with that induced by 
charge fluctuations. 7-121 

Note that the FS is already modified within the Hartree- 
Fock (HF) approximation, in which the self-energy does not 
include any fluctuation and is taken as = (-1/2) f G(k— 
?)(v c (?) + v s (q)). In the HF approximation, the FS is modi- 
fied simply by changing the effective hopping integrals f l; = 
tij — Yio- Vij( c i C/cr+h.c.), where t,j and Vy are the hopping in- 
tegral and Coulomb interaction between sites / and j, respec- 
tively. Thus, in the HF approximation the FS retains its round 
shape, whereas the FS in the FLEX approximation nests better 
because of its flat part, which appears only when fluctuations 
are fed back into 

To gain a deeper understanding of how the FS deforma- 
tion depends on charge fluctuations, we expand Eq. (6) with 
ft = //o + Sfi and kp = kpfi + SkpApfl, where //n and &f,o 
are the chemical potential and the Fermi wave vector for the 
noninteracting system, respectively, and Mf,o = &f,o/I&f,oI is 
a unit vector specifying a position on the FS. 7) A shift in the 
chemical potential due to Coulomb repulsion is then given by 
d/u / d 2 k 6(s k - ^ )Re £ s (&, 0)// d 2 k 6(e k - which cor- 
responds to an average of Re "E R (k, 0) on the noninteracting 
FS. Similarly, a shift of FS is given by 

fy-ReE s (A: Fi o,0) 



6k F 



IvroI 



(7) 



where Vf,o = dB/ c /dk\i [= i CF0 is the Fermi velocity. Here we 
used the fact that the direction of Vf,o is always directed op- 
posite to kp t o f° r the present hole-type FS. From Eq. (7), the 
deformation direction of FS (i.e., the sign of Skp) is deter- 
mined by the relative magnitude of ReX s (£, 0) on the FS 
of the noninteracting system (k = kp$). I n Fig- 3 (b), we 
show a contour plot of Re£ R (£, 0) for V = 2.43, where the 
FSs for the interacting and noninteracting systems are shown 
by the solid line and dashed line, respectively. The results 
show that the deformation of the FS is mainly induced by the 
relative increase of the exchange-correlation energy to holes 
(i.e., the decrease in Re I. R (k, 0) in terms of electrons) around 
k = (n/3, n/3) = Q/2, which corresponds to a hot spot suffer- 
ing strong electron scattering because of charge fluctuations. 

As a result of FS deformation, the quasiparticle-scattering 
anisotropy strengthens because the anisotropic self-energy is 
enhanced near CO. In Fig. 4, we show the quasiparticle scat- 
tering rate Im I. R (kp, 0) as a function of kp. The inset of Fig. 4 
shows three FS points: the -k x direction (A), +k y direction 
(B) and +k x direction (C). The results show that Im T. R (kp, 0) 
is largely enhanced in the hot region BC (the flat part of the 
FS) as the system approaches the CO transition, whereas it 
remains small in the cold region AB. Formation of hot and 
cold regions reflects anisotropy in the exchange-correlation 
interaction potential, which becomes strong near the CO tran- 
sition. 

This anisotropic quasi-particle scattering can be observed 
via the in-plane anisotropy of the electron transport such as 
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Fig. 4. Plot of -Im l. R (kp, 0) as a function of the position on the FS. Inset 
shows the shape of the FS for V = 2.4 in the region of k y > 0. Three 
representative positions on the FS are denoted by A, B, and C. 



(a) A(k,u) at V = 2 
4- — ■ 1 4 




7T j 2 



1.5 
1 

0.5 




Qo)A(k,u) at V = 2.4 

4i 



\ 




IT/2 



Fig. 5. Spectral weight A(k,a>) calculated for (a) V = 2 and (b) V = 2.4. 
The left (right) panels are contour plots of A(k, to) along the line k x = ky 
(k x = -k y ). Solid white lines are the band dispersions determined by oj = 
s k + Re I. R (k, to). 



electronic and thermal conductivities. Another experimen- 
tal probe suitable for observing both the FS deformation 
and anisotropic quasiparticle scattering is magnetoresistance. 
For example, out-of -plane conductivity cr :z under an in-plane 
magnetic field B may be calculated from the Boltzmann equa- 
tion in the relaxation-time approximation: 32-35 ' 

dk F 



0~ 7! oc t 



v F | 1 +(Q tF T iF ) 2 ' 



(8) 
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where t z (<K t) is the hopping integral between conduction 
planes, 1 /t> f = 2 |ImZ s (^F, 0)| is the quasiparticle lifetime, Vf 
is the Fermi velocity and = |vp x B|/|£fI is the cyclotron 
frequency. Eq. (8) indicates that the dependence of <x„ on the 
orientation of the in-plane magnetic field reflects the shape of 
the FS (through Vf and k-p) as well as the anisotropy in the 
quasiparticle lifetime. 36 ' 37 ' The calculation presented in this 
paper thus indicates that magnetoresistance in organic con- 
ductors is significantly affected by charge fluctuations near 
the CO transition. 

Finally, we discuss the spectral weight A(k, co) = 
-lmG R (k,u>)/7T, which can be measured directly by 
angle-resolved photoemission spectroscopy (ARPES) experi- 
ments. 39 ' We show the spectral weights for V — 2 and V = 2.4 
in Figs. 5 (a) and 5 (b), respectively. The spectral weight on 
the line k x — k y (in the hot region) is largely suppressed near 
the Fermi energy when the system approaches the CO transi- 
tion, while that on the line k x — —k y (in the cold region) near 
the Fermi energy remains unchanged. Note that, below the 
Fermi energy, a kink structure appears on the line k x = —k y . 
We attribute this high-energy kink structure to high-energy 
charge fluctuations. 40 ' 

We end this paper with a brief remark on CO phenom- 
ena in #-ET salts. In actual experiments on 9-ET salts, the 
horizontal-type CO, which is characterized by the wave vec- 
tor gh = (t/2, 7r/2) in our notation, is stabilized below the CO 
transition temperature. 41 ' The horizontal-type CO cannot be 
reproduced by the RPA (or mean-field) calculation on the ba- 
sis of the simple extended Hubbard model, 24 ' 25 ' and its origin 
has been discussed theoretically by considering the electron- 
phonon interaction. 42 ' Although the horizontal-type CO can- 
not be obtained in our calculation, we expect that our con- 
clusion on the FS deformation holds true also for this type of 
CO, because the rapid increase of V c {q) toward the CO tran- 
sition and the disagreement between the charge pattern and 
the nesting vector are essential to our conclusion, regardless 
of the detailed mechanism of the CO phenomena. 

In summary, we have studied the shape of the FS near 
the CO transition within the FLEX approximation. We have 
shown that, as the system approaches the CO transition, the 
FS is remarkably modified from that obtained from band cal- 
culations. This phenomenon is induced by strong charge fluc- 
tuations with a wave vector which is a large mismatch with 
the nesting vector of the FS for a noninteracting system. On 
the resultant FS, there appears a hot and cold region with dif- 
ferent quasiparticle lifetimes in each region. Our result may 
be observed experimentally via ARPES as well as by electron 
transport measurements such as those of in-plane magnetore- 
sistance anisotropy. Note that when the system is close to the 
CO transition, the FS may be sensitive to changes in temper- 
ature and/or pressure because strength of charge fluctuations 
changes rapidly. Our calculation can also be extended to more 
general CO models including those of inorganic systems. 
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